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Analysis  of  Algorithms  Predicting  Blood:Air  and  Tissue:Blood  Partition  Coefficients  from 
Solvent  Partition  Coefficients  for  use  in  Complex  Mixture  Physiologically  Based 
Pharmacokinetic/Pharmacodynamic  Modeling 


INTRODUCTION 

Jet  fuels  are  complex  petroleum  based  mixtures.  Analysis  of  JP-8  jet  fuel  by  mass  spectrometry 
yields  a  spectrum  of  approximately  300  identifiable  peaks  (Rodgers  et  al.,  1999).  U.S.  military 
personnel  have  been  and  continue  to  be  exposed  frequently  to  vapor,  aerosol  and  liquid  forms 
of  JP-8  since  it  became  the  standardized  military  fuel  (Olsen  et  al.,  1998;  Kobayashi  and 
Kikukawa,  2000).  To  make  the  situation  more  complicated,  military  exposures  to  JP-8  can 
occur  concurrently  with  exposures  to  high  levels  of  N,N-diethy!-m-toluamide  (the  insect  repellent 
DEET),  pyridostigmine  bromide  (anti-nerve  gas  agent)  and  permethrin  (pesticide),  which  may  be 
used  during  deployment.  Concurrent  military  occupational  exposures  have  been  implicated  in 
Gulf  War  Syndrome  (Peden-Adams  et  al.,  2001;  Riviere  et  al.,  2002). 

Low-level  JP-8  exposure  can  also  occur  through  environmental  pathways.  JP-8  spills  and 
leaking  underground  storage  tanks  contaminate  the  soil.  Of  the  1 .5  million  estimated  leaking 
underground  tanks  in  the  U.S.  in  1993,  40%  are  presumed  to  have  leaked  to  the  groundwater 
(Rodgers  et  al.,  1999).  Petroleum  fuels  are  contaminants  of  concern  at  an  estimated  50  to  60% 
of  all  Air  Force  cleanup  sites  (Miller,  1994).  Weathering  of  spilled  fuel  decreases  the  identifiable 
mass  spectrum  peaks  to  about  200,  still  a  very  complex  mixture  (Rodgers  et  al.,  1999).  Even 
more  complex,  people  environmentally  exposed  to  the  weathered  fuel  may  be  concurrently 
exposed  to  other  common  groundwater  contaminants  such  as  the  degreasing  solvents, 
trichloroethylene  and  perchloroethylene. 

Physiologically  based  pharmacokinetic  (PBPK)  modeling  of  mixtures  presents  a  challenge  for 
researchers.  A  preliminary  mixture  model  for  JP-8  in  the  rat  has  been  proposed  that 
interconnects  PBPK  models  of  individual  substances  via  competitive  metabolic  inhibition,  which 
produces  lower  overall  rates  of  elimination  with  increasing  chemical  complexity  (Robinson, 

2001 ).  Extrapolation  of  the  rat  model  to  humans  is  projected  for  the  future.  In  order  to 
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appropriately  test  and  validate  this  extrapolation,  blood:air  and  tissue:blood  partition  coefficients 
for  the  individual  components  are  required  to  describe  the  relative  distribution  of  each 
component  to  the  target  organ  (Gargas,  1991). 

The  partition  coefficient  (PC)  describes  the  distribution  of  a  chemical  between  two  media  at 
equilibrium.  A  tissue:b!ood  PC  is  necessary  for  each  compartment  in  the  model,  to  determine 
the  dose  distributed  by  the  rate  of  blood  flow  to  the  compartment,  which  represents  a  single 
tissue  or  lumped  (similar)  tissues.  A  blood:air  PC  is  necessary  for  volatile  chemicals  to 
determine  uptake  from  the  air  to  the  blood  during  inhalation  and  loss  of  chemical  to  the  air  from 
the  blood  during  exhalation  (Gargas,  1991). 

PCs  can  be  determined  experimentally  by  several  methods.  PCs  can  be  measured  in  vivo 
through  steady-state  exposure  (inhalation  or  continuous  iv  injection)  followed  by  rapid  analysis 
of  the  chemical  in  each  organ.  This  method  works  for  chemicals  that  are  not  metabolized  in  the 
body.  More  frequently,  PCs  are  determined  in  vitro.  Vial  equilibration  is  the  most  common  in 
vitro  method;  a  solvent  (oil,  octanol,  saline  or  water)  and  blood  or  tissue  (homogenized  or 
unhomogenized)  are  placed  in  the  bottom  of  a  small  vial.  The  vial  is  sealed  and  a  known 
amount  of  the  volatile  chemical  of  interest  is  injected.  The  vial  is  allowed  to  come  to  equilibrium; 
the  headspace  is  sampled  and  analyzed  on  a  gas  chromatograph.  This  method  allows 
deduction  of  the  concentration  in  the  medium  in  the  bottom  of  the  vial  without  having  to  analyze 
the  medium  directly  (Gargas,  1 991 ;  Gargas  et  at.,  1 989).  Adaptations  of  this  technique  exist  for 
nonvolatile  chemicals;  however,  many  of  these  adaptations  use  expensive  custom  radiolabeled 
chemicals  to  determine  the  concentration  of  the  chemical  at  equilibrium  between  the  tissue  and 
the  solvent  (oil  or  saline)  (Jepson  et  at.,  1994;  Murphy  et  al.,  1995).  The  cost  of  custom 
radiolabeled  chemicals  for  in  vitro  experiments  and  the  number  of  animals  needed  for  in  vivo  PC 
determination  are  reasons  why  nonvolatile  PCs  are  not  prevalent  in  the  literature  (Poulin  and 
Krishnan,  1995a).  Even  though  volatile  PCs  are  easier  to  determine,  tissue  and  blood  partition 
coefficients  are  not  available  for  all,  forcing  PBPK  modelers  to  look  for  alternatives  (Sweeney  et 
al.,  1996).  To  try  to  fill  this  void,  several  authors  have  developed  algorithms  as  ways  to  predict 
partition  coefficients  from  physical  properties. 

This  project  was  designed  to  examine  currently  available  methods  for  calculating  tissue  and 
blood  PCs  for  selected  JP-8  mixture  components  using  solvent  partition  coefficients  and  to 
compare  the  predicted  values  against  published  experimental  values.  In  this  context,  solvents 


2 


include  water  or  saline,  oil  and  octanol;  the  resulting  partition  coefficients  can  be  found  in  the 
literature  as  waterair,  saline:air,  oil:air  or  octanolrwater.  Additional  methods  of  determining  PCs 
include  quantitative  structure  activity  relationships  (QSAR)  calculations  which  require  extensive 
databases  of  physical  properties  (e.g.,  three-dimensional  structure,  octanoh.water  partition 
coefficients  (calculated  or  experimentally  derived),  normal  vapor  pressure)  for  hundreds  of 
chemicals  and  sophisticated  computer  programs  to  analyze  the  relationships  between  these 
parameters  and  the  dependent  variable  (blood :air  PC)  (Basak,  2002a,  2002b;  Gute  and  Basak, 
2001). 


METHOD 

Literature  Searches 

Literature  searches  were  performed  on  the  PubMed  (free  Medline  available  to  public)  and 
Science  Citation  Index  (Thompson  ISI,  Philadelphia,  PA)  databases  during  the  period  of  23  Jan 
through  12  Feb  2003.  The  searches  required  three  major  efforts.  First  solvent  PCs  were 
sought  using  the  word  “partition”,  the  solvent  types  (octanol,  water,  saline,  oil)  and  the  chemical 
names.  The  chemicals  selected  for  this  assessment  were  those  isolated  and  identified  in  a 
published  analysis  of  JP-8  by  gas  chromatography  (McDougal  et  al.,  2000).  The  14  most 
prevalent  chemicals  in  the  JP-8  sample  were  listed  along  with  their  concentration  (Table  1); 
isomers  (e.g.,  1,2,3-  vs.  1.2.4-trimethyl  benzene)  are  not  distinguishable  by  this  method. 
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Table  1.  Chemical  Composition  of  a  JP-8  Sample  as  Analyzed 
by  Gas  Chromatography  (McDougal  etal.,  2000) 


Chemical 

Undecane 

48.4 

Dodecane 

36.1 

Decane 

30.2 

Tridecane 

21.9 

Tetradecane 

14.6 

Methyl  naphthalenes 

9.9 

Trimethyl  benzenes 

9.7 

Nonane 

9.2 

Pentadecane 

8.4 

Dimethyl  naphthalenes 

6.3 

Xylenes 

4.8 

Naphthalene 

2.1 

Ethyl  benzene 

1.2 

Toluene 

0.5 

Second,  bloodrair,  tissue:air  and  tissue:blood  PCs  were  sought  using  the  word  “partition”,  blood 
or  tissue,  and  the  chemical  names.  Pharmacokinetic  models  for  all  the  chemicals  were  also 
searched,  as  any  PBPK  model  would  require  PCs  for  each  chemical.  The  third  literature  search 
effort  involved  locating  articles  on  calculating  PCs.  “Partition”  was  crossed  with  “predict”, 
“calculate”  and  “algorithm”.  Authors  found  to  be  prevalent  in  the  field  were  also  researched 
further  (e.g.,  Krishnan,  Basak).  Secondary  references  found  within  papers  were  sought  for  all 
three  literature  search  efforts. 


Experimental  PC  Tabulation 

Solvent  PCs  were  compiled  in  an  Excel  spreadsheet.  These  are  tabulated  in  Appendix  A,  Table 
A-1 .  Waterair  and  saline:air  partition  coefficients  were  combined  in  the  same  column. 

Saline:air  has  been  shown  to  equal  87  to  97%  of  the  waterair  PC  for  volatile  anesthetics.  Since 
this  difference  is  small,  introducing  little  error,  and  only  affects  the  most  volatile  compounds, 
saline:air  and  waterair  PCs  can  be  treated  as  the  same  (Meulenberg  and  Vijverberg,  2000). 

Experimental  tissue:air,  blood:air  and  tissue:blood  values  were  recorded  from  published  papers 
and  PBPK  models.  Results  were  published  as  either  tissue:air  or  tissuerblood  PCs.  Since 
PBPK  models  require  blood:air  and  tissue:blood  PCs,  not  tissue:air  PCs  (Gargas,  1991), 
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experimental  tissue:blood  PC  values  were  calculated  from  experimentally  derived  bloodrair  and 
tissue:air  PCs  by  the  following  equation  (Fiserova-Bergerova,  1983;  Poulin  and  Krishnan, 
1995a). 

Tissue  .Blood  Tissue 

- / - = -  (Equation  1) 

Air  Air  Blood 

PCs  were  separated  by  species  (rat  and  human)  and  averaged  by  tissue  and  chemical.  Some 
data  were  excluded,  as  discussed  in  the  Results  section.  The  final  values  used,  after 
exclusions  and  averaging,  are  listed  in  Tables  A-2  and  A-3. 

PC  Calculation 

Published  algorithms  to  predict  blood  and  tissue  partition  coefficients  using  solvent  PCs  were 
recorded  in  Tables  2  through  5  in  the  Results  section.  The  characteristics  of  the  equations  and 
the  authors’  correlation  coefficient  (R2)  were  recorded. 

Tissue  partition  coefficients  for  selected  JP-8  components  were  calculated  in  Excel 
spreadsheets.  The  number  of  chemicals  for  which  a  PC  was  calculated  for  any  given  algorithm 
was  limited  to  the  availability  of  the  solvent  PCs  necessary  to  perform  that  calculation;  if  a 
specific  solvent  PC  was  not  available  for  a  chemical  then  any  equation  using  that  solvent  PC 
could  not  be  worked  to  completion. 

Some  authors’  algorithms  directly  calculate  tissue:blood  PCs  (direct  equations)  while  the 
majority  result  in  tissue:air  values.  Equation  1  was  used  to  predict  the  tissuerblood  PC  from  a 
calculated  tissue:air  PC  value.  For  each  predicted  tissue:air  value,  tissue:blood  calculations 
were  performed  twice.  First,  an  experimentally  derived  blood:air  value  was  used  in  the 
denominator  of  Equation  1,  resulting  in  a  PCEb  value.  Second,  a  predicted  blood:air  value, 
using  the  appropriate  algorithm  from  the  same  author,  was  used  to  determine  a  tissue:blood 
PCpfi. 

Algorithm  results  generated  in  Excel  were  then  evaluated  in  SAS  (Statistical  Analysis  System, 
SAS  Institute,  Inc.,  Cary,  NC).  The  figures  in  this  report  were  prepared  in  SAS.  The  adequacy 
of  predictive  algorithms  was  evaluated  on  three  levels.  First,  the  performance  of  solvent  based 
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algorithms  as  a  whole  was  evaluated,  ranked  by  absolute  percent  error.  Twenty  percent  of  the 
measured  value  was  used  as  a  metric  for  overall  comparison.  Normal  biological  variation  in 
bloodrair  partition  coefficients  has  been  observed  at  20%  (Gearhart  et  al.,  1993).  Second,  each 
algorithm’s  average  performance  by  organ  was  compared  against  experimental  levels  using 
percent  error  plots  and  absolute  percent  error.  The  consistency  of  prediction  was  noted  (e.g., 
always  over-  or  under-predicts  compared  to  experimental  values)  and  the  adequacy  of 
prediction  of  alkane  PCs  (when  available)  was  considered.  Although  the  majority  of  JP-8  (Table 
1 )  consists  of  straight  chain  hydrocarbons,  the  most  data  are  available  for  the  aromatic 
constituents.  Third,  absolute  percent  error  was  used  to  compare  each  PCEb  and  PCPB  set  to 
determine  which  results  in  tissue:blood  values  closer  to  experimental  PCs.  The  minimum, 
median  and  maximum  relative  differences  between  PCEb  and  PCPB  values  were  found  using 
Equation  2. 


100x  PCEB  -PCPB 

Relative  Difference  =  — — - — -  (Equation  2) 

2 


The  percentage  of  occurrences  where  the  PCEB  value  had  a  smaller  absolute  percent  error  than 
the  PCPB  value  was  calculated. 


RESULTS 

Literature  Search  Results 

One  or  more  solvent  PCs  were  found  for  all  JP-8  components  sought.  Octanokwater  PCs  were 
found  for  all  compounds  in  Table  1;  handbooks  of  physical  chemistry  often  list  this  parameter 
(Yaws  etal.,  1999).  Octanokwater  PCs  are  most  often  measured  at  25°C.  Some  literature 
sources  did  not  state  the  temperature  at  which  this  parameter  was  measured,  but  the  values 
were  similar  to  or  the  same  as  values  measured  at  25°C  from  different  sources.  Therefore, 
when  more  than  one  octanokwater  PC  existed  for  a  chemical,  the  values  were  averaged.  One 
resource  measured  octanokwater  PC  at  37°C  (Poulin  and  Krishnan,  1995b);  this  value  was 
excluded  as  all  others  were  at  25°C  and  there  was  a  25°C  value  available  for  that  compound 
(toluene).  Temperature  purportedly  does  not  affect  octanokwater  partitioning  greatly  (Leo, 
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personal  communication,  2003).  The  toluene  log  octanohwater  PCs  were  2.65  and  2.73  at  37 
and  25°C,  respectively  (Gute  and  Basak,  2001;  Poulin  and  Krishnan,  1995b;  Yaws  etal.,  1999). 

The  other  solvent  PCs  (water  or  saline:air,  vegetable  or  olive  oikair)  were  harder  to  find.  This 
challenging  situation  has  been  noted  in  the  literature  (Poulin  and  Krishnan,  1995b).  These 
values  were  measured  at  37°C  in  all  cases. 

Empirical  tissue  and  blood  PCs  were  available  for  several  of  the  more  common  aromatics 
(toluene,  ethyl  benzene,  xylenes)  and  only  two  aliphatics  (decane  and  nonane).  The  best 
sources  for  these  parameters  were  papers  describing  methods  for  experimental  partition 
coefficient  determination  (e.g.,  Gargas  etal.,  1988,  1989;  Kumarathasan  etal.,  1998). 

Experimental  PC  Values 

When  multiple  tissue  and  blood  PC  values  were  found  for  a  chemical  in  a  tissue,  only  values 
from  primary,  experimentally  derived  sources  were  used  in  the  average.  All  primary  resources 
used  vial  equilibration  for  PC  determination  (Gargas  et  ai,  1989;  Kaneko  et  al.,  2000a; 
Kumarathasan  etal.,  1998;  Pierce  etal.,  1996;  Robinson,  2000).  Secondary  sources  did  not 
state  how  PCs  were  determined  (Fiserova-Bergerova,  1983;  Kaneko  et  al.,  2000b)  and  are 
always  subject  to  compounded  translation  errors. 

Other  values  that  were  excluded  from  the  averages  reported  in  Table  A-2  include  those  from 
older  rats.  Kumarathasan  etal.  (1998)  used  mature  (350  g)  male  Sprague-Dawley  rats  in 
several  of  their  experiments,  contrasting  their  results  with  those  of  250  g  young  adult  males  of 
the  same  strain.  PC  values  were  found  to  be  lower  in  the  younger  rats  and  more  comparable  to 
those  published  by  Gargas  etal.  (1989),  who  used  200  to  300  g  Fischer  344  male  rats 
(Kumarathasan  etal.,  1998).  It  is  expected  that  these  values  from  younger  rats  would  be  more 
useful  in  a  PBPK  model  as  young  adult  specimens  typically  are  used  in  kinetic  studies  for 
validation.  Age  increased  PCs  in  blood,  brain,  fat,  kidneys,  liver  and  muscle  by  an  average  of 
49,  160,  64,  240,  190  and  130%,  respectively,  for  m-  and  p-xylene.  The  authors  projected  that 
the  higher  PC  values  could  be  due  to  higher  fat  concentrations  or  altered  fat  composition  in  the 
older  animal  tissues  (Kumarathasan  et  al.,  1998).  Partition  coefficients  have  been  found  to  vary 
with  changes  in  hematocrit,  blood  lipids  (fasted  versus  postprandial  sampling)  and  organ  lipids 
(Fiserova-Bergerova,  1983). 


7 


Kumarathasan  et  al.  (1998)  reported  tissue:blood  PCs  instead  of  tissue:air  for  their  experiments 
with  250  g  rats.  For  m-  and  p-xylene,  the  authors  provided  a  b!ood:air  PC,  allowing  back- 
calculation  to  determine  tissue:air  PCs  compatible  for  comparison  with  values  from  other  studies 
(reverse  of  Equation  1).  There  were  no  biood:air  PCs  listed  for  ethyl  benzene  and  o-xylene. 
Young  adult  rat  blood:air  PCs  were  calculated  instead  from  mature  rat  PCs  reduced  by  49%,  the 
average  percent  difference  calculated  between  mature  and  young  adult  blood  PCs  for  m-  and  p- 
xylene.  For  o-xylene,  this  calculation  resulted  in  a  reduction  of  the  blood:air  PC  from  61.6  to 
41 .4,  a  similar  value  to  44.3,  reported  in  Gargas  et  al.  (1989).  A  comparable  value  for  ethyl 
benzene  was  not  found  but  the  calculation  reduced  the  PC  from  64.6  to  43.4  and  allowed  more 
credible  back-calculations  for  the  tissue:air  PCs. 

Fiserova-Bergerova  (1983)  reported  PCs  for  toluene  in  several  rat  tissues  (Table  A-2),  which 
are  consistently  lower  than  experimentally  derived  values  from  other  sources  (64,  30,  62  and 
35%  lower  than  experimental  values  for  blood,  fat,  kidney,  liver  and  muscle,  respectively). 
Because  experimental  details  were  not  available,  the  secondary  values  were  excluded  from 
averages  of  toluene  PCs  for  these  tissues.  Fiserova-Bergerova’s  secondary  values  for 
additional  tissues  not  found  in  other  studies  (brain,  heart,  kidney  and  lung)  remain  in  Table  A-2. 

Human  PCs  were  also  collected  and  listed  in  Table  A-3,  since  the  ultimate  goal  is  to  predict  the 
outcomes  of  exposure  to  mixtures  in  the  population  of  concern,  deployed  military  personnel 
(Robinson,  2001).  Sources  of  human  PCs  included  secondary  literature  references  (Fiserova- 
Bergerova,  1983;  Tardif  et  al.,  1995)  and  vial  equilibration  experiments  using  blood  from 
volunteers  or  fat  samples  from  patients  undergoing  surgery  (Gargas  et  al.,  1989;  Jarnberg  and 
Johanson,  1995;  Pierce  eta!.,  1996). 

Where  more  than  one  value  for  a  chemical  and  tissue  was  available,  only  results  from  similar 
experiments  and  techniques  were  averaged.  For  this  reason,  experiments  made  with 
perdeuterated  compounds  were  excluded  but  parallel  nondeuterated  PCs  were  used  in  the 
average.  Although  comparison  of  the  perdeuterated  chemical  PCs  with  the  nondeuterated  PCs 
revealed  that  they  were  not  statistically  different  for  all  chemicals  tested  (ethyl  benzene,  toluene, 
m-  and  o-xylene)  except  for  one  (p-xylene)  (Pierce  et  al.,  1996),  none  of  the  other  chemicals  for 
which  PCs  were  found  were  radiolabeled,  making  the  non-radiolabeled  numbers  more 
comparable. 
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A  fasted  value  for  toluene  was  excluded  in  the  average  blood:air  PC  in  humans.  During  fasting, 
blood  lipids  decrease,  resulting  in  a  lower  blood:air  PC  (Fiserova-Bergerova,  1983).  Since  the 
other  studies  did  not  specify  the  food  status  of  subjects,  it  is  assumed  that  blood  was  taken  in 
various  fasted  or  postprandial  states  depending  on  the  individual.  The  one  to  two  hour 
postprandial  value  (Fiserova-Bergerova,  1983)  was  comparable  to  values  from  different  sources 
and  was  kept  in  the  final  average  for  this  chemical. 

Back  calculations  to  derive  tissue:air  PCs  were  performed  using  the  reverse  of  Equation  1 . 
Tissue:blood  PC  values  were  multiplied  by  the  bloodrair  PC  value  reported  by  the  same  authors 
in  the  same  study  (e.g.,  Kumarathasan  et  al.  (1998)  fat:air  value  for  250  g  rats  was  divided  by 
the  blood:air  value  for  250  g  rats).  However,  blood:air  PCs  were  not  always  reported  in  each 
study.  For  example,  Fiserova-Bergerova  (1983)  did  not  report  a  blood:air  PC  for  toluene  in  rats, 
so  the  tissue:blood  values  were  multiplied  by  the  experimental  blood:air  value  reported  by 
Gargas  et  al.  (1989).  Calculations  to  determine  tissue:blood  PCs  from  tissue:air  and  blood:air 
PCs  were  completed  using  Equation  1  with  values  from  the  same  study  whenever  possible. 
These  calculated  values  are  listed  in  Tables  A-2  and  A-3. 

Calculated  PC  Results 

Several  solvent  PC  based  equations  were  found  to  predict  tissue  and  blood  PCs.  These 
equations  are  listed  in  Tables  2  through  5.  The  equations  were  categorized  as  rat  or  human 
and  according  to  the  result  of  the  equation  (tissueiair  or  blood:air  versus  tissue:blood).  Three 
types  of  equations  were  found.  Empirical  equations  are  derived  from  linear  regressions  with 
experimental  data.  These  equations  are  assumed  to  provide  best  estimated  values  for 
chemicals  similar  to  the  type  of  compounds  used  to  derive  the  equation.  For  example,  low 
molecular  weight  volatile  PCs  are  expected  to  be  predicted  better  by  Equation  5  than  Equation  6 
(see  Table  2)  as  PCs  from  this  type  of  compound  were  used  by  Gargas  et  al.  (1989)  to  develop 
Equation  5  (Payne  and  Kenny,  2002).  The  authors’  R2  values  were  included  in  the  table,  when 
available,  to  show  the  strength  of  relationship  between  the  experimental  data  used  to  derive  the 
equation  and  the  equation  predicted  values. 

Physiologically  based  equations  take  into  account  the  water  and  lipid  make-up  of  the  tissue  in 
question.  Poulin  and  Krishnan  (1996a)  (e.g.,  Table  2,  Equation  3)  describe  blood  and  other 
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tissues  as  being  made  up  of  tissue  water,  neutral  lipids  (non-polar  lipids  including  triglycerides 
and  cholesterol)  and  phospholipids  (lipids  with  a  phosphoric  acid  ester  on  the  glycerol  molecule 
including  phosphatidyl  choline).  Phospholipids  are  both  lipid  and  water  soluble.  Physiologically 
based  equations  do  not  contain  any  additional  empirical  components  and  are  not  fitted  through 
regression. 

Hybrid  equations  combine  elements  of  physiologically  based  and  empirical  equations  (e.g.. 
Table  2,  Equation  4).  The  authors  of  hybrid  equations  attempt  to  account  for  lipophilic  and 
hydrophilic  components  of  each  tissue,  and  then  add  an  empirical  component  to  adjust  for 
systematic  errors  in  the  fit  of  the  equation  or  to  represent  protein  binding  in  the  blood  or  tissue 
(Meulenberg  and  Vijverberg,  2000).  Like  strict  empirical  equations,  the  predicted  PCs  from 
hybrid  equations  are  expected  to  be  more  accurate,  for  example,  for  volatiles  if  they  were 
derived  from  volatiles  and,  therefore,  are  not  expected  to  work  as  well  for  non-volatiles 
(DeJongh  et  al.,  1997).  The  authors’  R2  values  were  included  in  the  table,  when  available. 
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Table  2  Rat  Blood:Air  and  TissuerAir  Algorithms 
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As  there  were  few  sources  for  solvent  PC  based  equations  for  direct  human  tissue:blood  PCs 
(Table  5),  fatrblood  PC  based  equations  by  Fiserova-Bergerova  and  Diaz  (1986)  were  included 
(e.g.,  Table  5,  Equation  44).  The  small  number  of  chemicals  with  experimental  fat:blood  PCs  in 
humans  limited  the  number  of  predicted  PCs  calculable  for  this  set  of  equations. 

Overall  Performance  of  Algorithms 

Complete  results  for  all  algorithms,  tissues  and  chemicals  are  listed  in  Appendix  B.  The  percent 
errors  for  each  blood:air  and  PCEb  tissue:blood  calculation  for  rats  and  humans  are  listed  below 
in  Tables  6  and  7,  respectively. 
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Table  6:  Percent  Error  from  Rat  Blood:Air  and  Tissue:Blood  (Direct  and 

PCEb  Equations)  Algorithms 


Percent  Error 

Tissue 

Chemical 

D 

G8 

G9 

M 

P5 

P6 

Decane 

-83 

-13 

356 

-95 

■■ 

-16 

34 

-48 

-87 

Blood:Air 

Toluene 

31 

72 

-52 

-79 

Xylene,  m- 

-8 

43 

-46 

-83 

Xylene,  o- 

8 

60 

-44 

-80 

riii  ■l«i 

-8 

44 

-45 

-85 

■M 

480 

Toluene 

-23 

Brain:Blood 

Xylene,  m- 

294 

Xylene,  o- 

186 

Xylene,  p- 

306 

118 

120 

68 

41 

-96 

28 

Nonane 

-28 

-99 

Fat:Blood 

Toluene 

-10 

-5 

-17 

-36 

-97 

-8 

Xylene,  m- 

71 

66 

28 

8 

-97 

23 

Xylene,  o- 

34 

42 

9 

-7 

-97 

18 

70 

80 

39 

16 

-97 

11 

liUf.l, !«.!.%■ 

Toluene 

234 

104 

32 

132 

93 

-8 

93 

Nonane 

415 

-15 

Liver:  Blood 

Toluene 

-59 

-58 

-41 

-62 

-76 

-44 

Xylene,  m- 

37 

-9 

58 

28 

-40 

63 

2 

-22 

28 

2 

-54 

43 

Xylene,  p- 

26 

-12 

54 

25 

-43 

34 

ima 

-19 

-41 

-6 

-11 

-3 

-33 

Nonane 

5 

-30 

Muscle:Blood 

Toluene 

-54 

-45 

-32 

-55 

-37 

00 

1 

Xylene,  m- 

-14 

-34 

3 

-5 

1 

-9 

Xylene,  o- 

-30 

-39 

-9 

-19 

-16 

-15 

Xylene,  p- 

-24 

-40 

-5 

-11 

-8 

-29 

Notes:  D  =  DeJongh  etal.  (1997)  (H),  G8  =  Gargas  et  al.  (1988)  (E),  G9  =  Gargas  etal.  (1989)  (E),  M  = 
Meulenberg  and  Vijverberg  (2000)  (H),  P5  =  Poulin  and  Krishnan  (1995b)  (P),  P6  =  Poulin  and  Krishnan 
(1996a)  (P),  where  (E)  =  empirical,  (H)  =  hybrid  and  (P)  =  physiological 
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Table  7:  Percent  Error  from  Human  Blood:Air  and  Tissue:Blood  (Direct  and  PCEB 

Equations)  Algorithms 


Percent  Error 

Tissue 

Chemical 

A 

B 

D 

F 

G9 

M 

P6 

T 

Blood  Air 

Ethyl  benzene 

-47 

27 

2 

-61 

-72 

Toluene 

-21 

26 

-34 

-59 

-49 

Trimethyl  benzene,  1,2,3- 

-56 

17 

22 

-21 

-78 

Trimethyl  benzene,  1,2,4- 

-66 

6 

27 

-32 

-85 

Trimethyl  benzene,  1,3,5- 

-62 

31 

68 

-67 

-83 

Xylene,  m- 

-52 

9 

-16 

-60 

-74 

Xylene,  o- 

-41 

22 

-12 

-54 

-66 

Xylene,  p- 

-65 

-18 

-36 

-74 

-81 

Brain:Blood 

Toluene 

-46 

-19 

-10 

-8 

-55 

-34 

Fat:Blood 

Ethyl  benzene 

41 

-61 

70 

-4 

8 

7 

Toluene 

-22 

-75 

28 

-54 

-18 

-31 

Xylene,  m- 

22 

-59 

77 

-18 

17 

-6 

Xylene,  o- 

9 

-66 

50 

-28 

13 

-13 

Xylene,  p- 

7 

-53 

104 

-27 

-13 

-18 

Heart:  Blood 

Toluene 

21 

Kidney:  Blood 

Toluene 

-49 

-13 

33 

26 

-47 

-38 

Liver:  Blood 

Toluene 

-72 

-41 

8 

-21 

-22 

-48 

Xylene,  m- 

-67 

34 

122 

4 

40 

-13 

Lung:  Blood 

Toluene 

-85 

-50 

-31 

-63 

Muscle:Biood 

Toluene 

-68 

-39 

-3 

-17 

-66 

-13 

Note:  A  =  Abraham  et  al.  (1985)  (E),  B  =  Balaz  and  Lukacova  (1999)  (H),  D  =  DeJongh  et  at.  (1997)  (H), 
F  =  Fiserova-Bergerova  and  Diaz  (1986)  (E),  G9  =  Gargas  et  al.  (1989)  (E),  M  =  Meulenberg  and 
Vijverberg  (2000)  (H),  P6  =  Poulin  and  Krishnan  (1996a)  (P),  T  =  Tichy  (1991)  (E) 


The  values  in  Tables  6  and  7  vary  widely  (-99  to  480%)  and  do  not  give  a  good  initial  impression 
of  algorithms  used  to  calculate  PCs.  Table  8  provides  perspective  for  evaluating  these 
equations.  Use  of  Table  8  is  illustrated  in  the  following  example:  for  all  blood:air  and 
tissue:blood  (direct  and  PCEB  equations)  algorithms,  29%  of  the  calculated  values  were  within 
±20%  of  the  experimental  value.  Although  Table  8  lists  additional  important  percentages, 
attention  is  called  to  the  ±20%  minimum  standard  set  in  the  Methods  section  due  to  normal 
levels  of  variability  in  biological  data  (Gearhart  etal.,  1993). 
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Table  8:  Percent  of  Blood:Air  and  Tissue:Blood  (Direct  and  PCEb  Equations)  Algorithms 

within  ±  Absolute  Percent  Error  Levels 


1 

Number  of 
Calculations 

Absolute  Percenl 

Error 

10 

20 

30 

40 

50 

90 

All 

231 

15 

29 

42 

54 

65 

90 

93 

19 

30 

42 

56 

70 

98 

■EBH 

83 

12 

28 

42 

53 

63 

86 

■-j  wmmmw 

55 

13 

31 

40 

51 

60 

86 

The  type  of  equation  did  not  affect  the  percentage  of  the  calculations  with  values  within  20%  of 
experimental  PCs;  all  types  had  nearly  30%  of  outcomes  in  this  category.  The  type  or  basis  of 
the  algorithm  appears  to  have  little  effect  on  accuracy  of  the  outcome  and  no  type  of  equation 
can  be  considered  as  better  than  another. 

Performance  by  Tissue  and/or  Chemical 

The  algorithms  were  next  evaluated  on  the  basis  of  prediction  adequacy  for  each  tissue  type. 

Blood: Air  -  Figures  1  and  2  depict  the  spread  of  predictions  around  the  average  experimental 
value  for  rat  and  human  blood:air  PCs,  respectively.  Each  figure  lists  the  average  for  each 
algorithm  source  at  the  right  side.  As  there  was  a  predicted  value  for  decane  rat  blood:air  PC, 
and  as  decane  was  the  only  alkane  for  which  experimental  values  existed  to  facilitate 
comparison,  two  sets  of  averages  are  listed  for  Figure  1,  with  and  without  decane. 
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Figure  1 :  Percent  Error  for  Each  Estimate  of  Rat  Blood:  Air  PCs 
Relative  to  Experimental  Values 
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80- 

G9 

with  Decane 

60- 

G9 

G8  =  26 

k— 

o 

40- 

G9 

G8 

G9 

G9 

G9  =  44 

M  =98 

L_ 

LU 

20- 

G8 

P6  =  85 

without  Decane 

c 

<D 

O 

i_ 

-20- 

G9 

G8 

G8 

G8 

G8  =  14 

G9  =  51 

0) 

a. 

-40- 

M 

M 

M 

M 

M 

M  =47 

P6  =83 

-60- 

-80- 

G8 

P6 

P6 

P6 

P6 

P6 

”1” 

O 

o 

P6 

1 — i - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 — 1 

Decane  Ethyl  Toluene  Trimethyl  Trimethyl  Trimethyl  Xylene,  Xylene,  Xylene, 
benzene  benzene,  benzene,  benzene,  m-  o-  p- 

1,2,3-  1,2,4-  1,3,5- 


Notes:  Circled  letters  =  out  of  range  (greater  than  100%  error),  G8  =  Gargas  et  al.  (1988)  (E),  G9  = 
Gargas  et  al.  (1989)  (E),  M  =  Meulenberg  and  Vijverberg  (2000)  (H),  P6  =  Poulin  and  Krishnan  (1996a) 
(P),  where  (E)  =  empirical,  (H)  =  hybrid  and  (P)  =  physiological 


Figure  2:  Percent  Error  for  Each  Estimate  of  Human  Blood:Air  PC 
Relative  to  Experimental  Values 


Average  |%Error| 

A  =51 
G9  =  20 
M  =27 
P6  =54 
T  =74 


Notes:  A=  Abraham  et  al.  (1985)  (E),  G9  =  Gargas  et  al.  (1989)  (E),  M  =  Meulenberg  and  Vijverberg 
(2000)  (H),  P6  =  Poulin  and  Krishnan  (1996a)  (P),  T  =  Tichy  (1991)  (E),  where  (E)  =  empirical,  (H)  = 
hybrid  and  (P)  =  physiological 


Predicted  rat  blood:air  PCs  were  best  described  by  Gargas  etal.  (1988)  (G8,  Equation  6).  The 
average  absolute  percent  error  for  this  source  was  26%  with  decane  (14%  without).  In  addition, 
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values  from  this  source  had  negative  percent  errors  for  some  chemicals  and  positive  for  others, 
indicating  that,  for  any  new  chemical  similar  to  the  ones  used  in  this  report,  the  result  of  this 
equation  should  fall  near  the  experimental  value  and  not  be  heavily  skewed  like  some  equations 
(P  and  M,  Equations  3  and  4),  resulting  in  consistent  severe  underprediction  (Poulin  and 
Krishnan,  1996a  and  Meulenberg  and  Vijverberg,  2000,  respectively). 

Average  absolute  percent  error  values  decreased  for  three  of  the  four  sources  of  rat  bloodrair 
PC  algorithms  when  the  predicted  decane  value  was  removed.  Only  the  Gargas  etal.  (1989) 
equation  (G9,  Equation  5)  adequately  predicted  the  bloodrair  PC  for  this  alkane.  Although  the 
Gargas  etal.  (1989)  equation  performed  better  when  alkanes  were  included,  its  overall 
averages  were  higher  (|44%|  with  decane)  than  Gargas  et  al.  (1988)  (G8,  Equation  6,  |26%|  with 
decane)  and  so  would  not  be  the  best  algorithm  for  this  chemical  set.  Given  the  highly  alkane 
composition  of  JP-8  (McDougal  et  al.,  2000),  the  algorithms  evaluated  here  may  not  be 
appropriate  for  this  task. 

In  comparing  predicted  rat  and  human  bloodrair  PC  distributions  around  the  experimental  values 
(Figures  1  and  2),  the  spread  of  percent  error  values  is  wider  for  predicted  rat  PCs  while  the 
spread  for  human  percent  errors  is  somewhat  smaller.  The  differing  success  in  predicting  rat 
and  human  values  is  apparent  when  comparing  equations  from  the  same  source  (e.g.,  Gargas 
et  al.  (1989),  G9,  Equations  5  and  23).  The  greater  difficulty  in  calculating  rat  blood:air  PCs 
may  be  related  to  increased  binding  to  proteins  occurring  in  rat  blood  as  compared  to  humans 
(Poulin  and  Krishnan,  1996b).  Protein  binding  may  also  be  a  factor  in  why  the  physiological 
equations  of  Poulin  and  Krishnan  (1996a)  (P6,  Equations  3  and  21)  performed  more  poorly  in 
predicting  PCs  for  this  medium  (as  compared  to  the  empirical  equations),  as  their  approach 
does  not  incorporate  a  binding  or  correction  term. 

The  human  blood:air  experimental  PC  values  were  best  predicted  by  Gargas  etal.  (1989)  (G9, 
Equation  23),  with  an  average  absolute  percent  error  of  20%.  The  Meulenberg  and  Vijverberg 
(2000)  equation  (M,  Equation  22)  also  did  reasonably  well  with  an  average  absolute  value  of 
27%.  Since  experimental  values  for  decane  human  blood:air  PCs  were  not  found,  a 
comparison  of  predictions  with  and  without  an  alkane  was  not  possible. 

Brain:Blood  -  Only  one  source,  Meulenberg  and  Vijverberg  (2000)  (Equation  7),  was  found  to 
provide  a  solvent-based  brain:air  equation  with  which  to  calculate  a  rat  brainrblood  PCEb- 
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Therefore,  no  figure  for  comparison  can  be  constructed.  As  seen  in  Table  6,  percent  error  for 
individual  estimates  as  compared  to  the  experimental  value  varied  widely  from  -23  to  480%;  it  is 
unlikely  that  this  equation  would  be  considered  useful  for  modeling  purposes  for  this  set  of 
chemicals.  Conversely,  only  a  single  experimental  value  was  found  for  blood:brain  PCs  in 
humans.  Six  different  sources  were  located  with  algorithms  for  this  parameter.  With  only  a 
single  predicted  value  each  to  compare  with  an  experimental  value  (Table  7),  no  meaningful 
comparison  of  the  algorithms  can  be  discussed. 

Fat:Blood~  Figures  3  and  4  show  the  individual  predictions  around  the  experimental  value  for 
rat  and  human  fat:blood  PCs,  respectively.  These  PCs  were  derived  either  from  direct 
equations  (Tables  4  and  5)  or  calculated  from  predicted  tissue:air  and  experimental  blood:air 
values  (PCEb,  Tables  2  and  3).  Averages  for  rat  PCs  again  were  listed  with  and  without  the 
available  alkane,  nonane.  Human  experimental  values  were  not  available  for  this  chemical. 


Figure  3:  Percent  Error  (Relative  to  Experimental  Values)  for  Each 
Estimate  of  Rat  Fat:Blood  PC  (Direct  and  PCEB  Equations) 


Average  |%Error| 

with  Nonane 
D  =55 
G8  =  63 
G9  =  32 
M  =21 
P5  =97 
P6  =  18 


Notes:  Circled  letters  =  out  of  range  (greater  than  100%  error);  D  =  DeJongh  etal.  (1997)  (H);  G8  = 
Gargas  etal.  (1988)  (E);  G9  =  Gargas  etal.  (1989)  (E);  M  =  Meulenberg  and  Vijverberg  (2000)  (H);  P5  = 
Poulin  and  Krishnan  (1995h)  (P);  PR  =  Pnulin  and  Krishnan  (1996a)  (P),  where  (E)  =  empirical.  (H)  = 
hybrid,  (P)  =  physiological 
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Figure  4:  Percent  Error  (Relative  to  Experimental  Values)  for  Each  Estimate 
of  Human  Fat:Blood  PC  (Direct  and  PCEB  Equations) 


Average  |%Error| 

A  =20 
B  =63 
D  =66 
M  =26 
P6  =  14 
T  =15 


Notes:  Circled  letters  =  out  of  range  (greater  than  100%  error),  A  =  Abraham  et  al.  (1985)  (E);  B  =  Balaz 
and  Lukacova  (1999)  (H);  D  =  DeJongh  et  al.  (1997)  (H);  M  =  Meulenberg  and  Vijverberg  (2000)  (H);  P6  = 
Poulin  and  Krishnan  (1996a)  (P);  T  =  Tichy  (1991)  (E),  where  (E)  =  empirical,  (H)  =  hybrid,  (P)  = 
physiological 


For  both  rat  and  human  fat:blood  PCs,  Poulin  and  Krishnan  (1996a)  (P6,  Equations  8  and  29) 
performed  best.  Average  absolute  percent  error  in  rat  fat:blood  PC  predictions  was  18%;  since 
these  are  physiological  equations,  the  tissue  components  appear  to  be  most  important  in 
predicting  the  partitioning  of  chemicals  in  fat.  In  human  predictions,  Poulin  and  Krishnan 
(1996a)  values  averaged  only  14%  absolute  percent  error.  Meulenberg  and  Vijverberg  (2000) 
values  (M,  Equations  9  and  30)  also  performed  well  in  both  species,  with  average  absolute 
percent  errors  of  21  and  26%  for  rat  and  human  fat,  respectively.  Meulenberg  and  Vijverberg 
(2000)  used  hybrid  equations  that  incorporated  terms  for  tissue  makeup.  In  humans,  the  Tichy 
(1991)  equation  (T,  Equation  27)  and  the  Abraham  et  al.  (1985)  (A,  Equation  32)  also  performed 
well  at  15%  and  20%  absolute  percent  error,  respectively.  Abraham  et  al.  and  Tichy  both 
employed  empirical  equations  generated  from  linear  regressions. 

Appropriate  solvent  PCs  were  not  available  for  decane  to  use  all  of  the  rat  algorithms;  only  the 
two  equations  using  octanokwater  partitioning  could  be  calculated  for  this  compound. 

Therefore,  the  efficacy  of  these  equations  for  alkane  chemicals  cannot  be  ascertained.  No 
human  experimental  PCs  for  decane  were  found. 
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Heart: Blood  -  No  heart:blood  PC  algorithms  were  found  for  rats.  A  single  direct  equation  for 
heartblood  PC  prediction  was  located  for  humans  (Balaz  and  Lukacova,  1999,  Equation  55). 

An  experimental  value  for  heart:blood  PC  was  found  for  one  chemical,  toluene.  Therefore  a 
single  prediction  could  be  calculated,  resulting  in  a  21%  error  as  compared  to  the  experimental 
number.  No  further  comparisons  can  be  made. 

Kidney: Blood  -  A  single  source,  Meulenberg  and  Vijverberg  (2000)  (Equation  12),  was  found  to 
provide  a  solvent-based  kidneyrair  equation  with  which  to  calculate  a  rat  kidney:blood  PCEB. 
Compared  with  a  single  experimental  value  for  toluene,  this  prediction  was  234%  different 
(Table  6).  No  further  comparison  can  be  constructed.  Six  sources  were  located  providing 
algorithms  for  human  kidneyrair  PCs.  Again,  only  a  single  experimental  value  was  found  for 
blood:brain  PCs  in  humans.  With  only  a  single  predicted  value  each  to  compare  with  an 
experimental  value  (Table  7),  no  meaningful  comparison  of  the  algorithms  can  be  discussed. 

LiverBlood  -  Figure  5  depicts  the  individual  predictions  around  the  experimental  value  for  rat 
liver:blood  PCs.  These  PCs  were  either  direct  or  PCEb  calculated. 
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Figure  5:  Percent  Error  (Relative  to  Experimental  Values)  for  Each  Estimate 
of  Rat  Liver:Blood  PC  (Direct  and  PCEB  Equations) 


Average  |%Error| 

with  Nonane 
D  =107 
G8  =  26 
G9  =  62 
M  =  42 
P5  =  39 
P6  =  55 


Notes:  Circled  letters  =  out  of  range  (greater  than  100%  error),  D  =  DeJongh  et  al.  (1997)  (H),  G8  = 
Gargas  etal.  (1988)  (E),  G9  =  Gargas  et  al.  (1989)  (E),  M  =  Meulenberg  and  Vijverberg  (2000)  (H),  P5  = 
Poulin  and  Krishnan  (1995b)  (P),  P6  =  Poulin  and  Krishnan  (1996a)  (P),  where  (E)  =  empirical,  (H)  = 
hybrid,  (P)  =  physiological 


For  rat  liver:blood  PC  predictions,  Gargas  et  al.  (1988)  (G8,  Equation  16)  had  the  lowest 
average  absolute  percent  error  (26%)  across  the  chemicals.  Again,  only  the  two  equations 
using  octanol:water  partitioning  could  be  calculated  for  this  compound.  Therefore,  the  efficacy 
of  these  equations  for  alkane  chemicals  can  not  be  ascertained.  Only  two  human  experimental 
values  for  liverblood  PC  were  found,  so  a  figure  was  not  generated.  With  just  two  compounds 
for  comparison  (toluene  and  m-xylene),  averages  were  not  calculated.  A  study  of  Table  7 
indicates  that  Fiserova-Bergerova  and  Diaz  (1986,  Equation  60)  would  have  the  lowest  absolute 
percent  error;  however  the  experimental  data  are  too  limited  to  justify  this  conclusion. 

Lung:Blood  -  Lung:blood  PC  algorithms  were  found  only  for  humans.  An  experimental  PC  was 
located  for  toluene.  The  four  algorithms  ranged  from  -31  to  -85%  error  (Table7).  Further 
conclusions  cannot  be  made. 

Muscle:Blood  -  Figure  6  shows  the  spread  of  experimental  values  for  rat  muscle:blood  PCs. 
Again,  these  PCs  were  either  direct  or  PCEB  calculated. 
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Figure  6:  Percent  Error  (Relative  to  Experimental  Values)  for  Each  Estimate 
of  Rat  MusclerBlood  PC  (Direct  and  PCEb  Equations) 


Average  |%Error| 

with  Nonane 
D  =24 
G8  =  40 
G9=  11 
M  =20 
P5  =  16 
P6  =27 


Notes:  D  =  DeJongh  et  al.  (1997)  (H),  G8  =  Gargas  et  al.  (1988)  (E),  G9  =  Gargas  et  at.  (1989)  (E),  M  = 
Meulenberg  and  Vijverberg  (2000)  (H),  P5  =  Poulin  and  Krishnan  (1995b)  (P),  P6  =  Poulin  and  Krishnan 
(1996a)  (P),  where  (E)  =  empirical,  (H)  =  hybrid,  (P)  =  physiological 


Overall,  rat  muscle:blood  PC  tended  to  be  underpredicted  by  algorithms.  Several  sources  did 
adequate  jobs  of  predicting  this  PC,  judging  by  the  average  absolute  percent  errors,  including 
Gargas  etal.  (1989)  (G9,  Equation  19)  at  11%,  Poulin  and  Krishnan  (1995b)  (P5,  Equation  50) 
at  16%  and  Meulenberg  and  Vijverberg  (2000)  (M,  Equation  18)  at  20%. 

Human  experimental  values  for  muscle:blood  PC  were  limited  to  one,  toluene.  Six  sources 
provided  algorithms  to  calculate  this  PC.  Percent  errors  ranged  from  -3  to  -68%  (Table  7). 

A  synopsis  of  conclusions  that  may  be  drawn  following  the  tissue  specific  assessment  of  PC 
prediction  algorithms  can  be  found  in  Table  9  below.  There  were  few  tissues  (four  rat  and  two 
human)  for  which  both  experimental  PC  values  and  equations  could  be  located  to  allow  for 
meaningful  comparison  of  the  predictive  algorithms. 
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Table  9.  Synopsis:  Best  Performance  of  Algorithms  by  Species  and  Organ 

for  Solvent  PC  Based  Equations 


Rat 

Human 

PC 

Best  Equation  for  Chemical  Set 

Based  on  Average  |%  Error| 

Best  Equation  for  Chemical  Set 

Based  on  Average  |%  Error| 

Blood:Air 

Gargas  et  al.  (1988) 

Equation  6 

|26%|  including  alkanes 
|14%|  excluding  alkanes  • 

Gargas  ef  al.  (1989) 

Equation  23 
|20%| 

Meulenberg  and  Vijverberg  (2000) 

Equation  22 
|27%| 

Brain:Blood 

insufficient  equations  for  comparison 

insufficient  experimental  data  for 
comparison 

Fat:  Blood 

Poulin  and  Krishnan  (1996a) 

Equation  8 

1 1 8%| 

Meulenberg  and  Vijverberg  (2000) 

Equation  9 
|21  %| 

Poulin  and  Krishnan  (1996a) 

Equation  29 
|14%| 

Tichy  (1991) 

Equation  27 
|1 5%| 

Abraham  ef  al.  (1985) 

Equation  32 
|20%| 

Meulenberg  and  Vijverberg  (2000) 

Equation  30 
|26%| 

Heart:Blood 

no  equation  found 

insufficient  equations  and 
experimental  data  for  comparison 

Kidney:  Blood 

insufficient  equations  for  comparison 

insufficient  experimental  data  for 
comparison 

Liver:  Blood 

Gargas  (1988) 

Equation  16 
|26%| 

insufficient  experimental  data  for 
comparison 

Lung:Blood 

no  equation  found 

insufficient  experimental  data  for 
comparison 

Muscle  Blood 

Gargas  ef  al.  (1989) 

Equation  19 
|11%| 

Poulin  and  Krishnan  (1995b) 

Equation  50 
|16%| 

Meulenberg  and  Vijverberg  (2000) 

Equation  18 
|20%| 

insufficient  experimental  data  for 
comparison 

PCs  for  Toluene  -  To  allow  for  a  different  perspective,  a  profile  of  predictions  for  a  single 
chemical  was  considered  for  both  rat  and  human  tissues  (Figures  7  and  8,  respectively). 


Toluene  was  chosen  as  it  had  the  most  complete  list  of  experimental  PC  values.  Overall 
performance  was  evaluated  using  the  average  absolute  percent  error  for  each  source  of 
algorithms;  these  percents  are  listed  to  the  right  of  the  figures. 


Figure  7:  Percent  Error  (Relative  to  Experimental  Values)  for  Each  Estimate  of  Rat 
Tissue:Blood  PC  (Direct  and  PCEb  Equations)  for  Toluene 


Rat  Toluene 


Average  |%Error| 

D  =41 
G8  =  36 
G9  =  30 
M  =82 
P5  =  70 
P6  =  33 


Notes:  Circled  letters  =  out  of  range  (greater  than  100%  error),  D  =  DeJongh  et  al.  (1997)  (H),  G8  = 
Gargas  et  al.  (1988)  (E),  G9  =  Gargas  et  al.  (1989)  (E),  M  =  Meulenberg  and  Vijverberg  (2000)  (H),  P5  = 
Poulin  and  Krishnan  (1995b)  (P),  P6  =  Poulin  and  Krishnan  (1996a)  (P),  where  (E)  =  empirical,  (H)  = 
hybrid,  (P)  =  physiological 
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Figure  8:  Percent  Error  (Relative  to  Experimental  Values)  for  Each  Estimate  of  Human 
Tissue:Blood  PC  (Direct  and  PCEB  Equations)  for  Toluene 


Human  Toluene 


Brain  Fat  Kidney  Liver  Muscle  Heart  Lung 


Tissue 


Notes:  A  =  Abraham  et  al.  (1985)  (E),  B  =  Balaz  and  Lukacova  (1999)  (H),  D  =  DeJongh  et  ah  (1997)  (H), 
F  =  Fiserova-Bergerova  and  Diaz  (1986)  (E),  M  =  Meulenberg  and  Vijverberg  (2000)  (H),  P6  =  Poulin  and 
Krishnan  (1996a)  (P),  T  =  Tichy  (1991)  (E),  where  (E)  =  empirical,  (H)  =  hybrid,  (P)  =  physiological 


A  different  picture  is  formed  when  looking  at  a  single  chemical  rather  than  by  tissue  (Table  9). 
Algorithms  from  Gargas  et  al.  (1989)  performed  “best”  (lowest  average  absolute  percent  error, 
30%)  across  rat  tissues,  while  DeJongh  et  al.  (1997)  was  “best”  for  human  tissues  (16%);  Poulin 
and  Krishnan  (1 996a)  was  the  second  “best  source”  for  both  species.  In  tissue  by  tissue 
comparisons  above,  Gargas  et  al.  (1989)  was  found  to  be  the  best  predictor  for  rat  muscle:blood 
PCs  but  was  not  otherwise  recommended  in  Table  9  for  calculation  of  rat  PC  values.  DeJongh 
et  al.  (1997)  was  not  recommended  in  Table  9  for  any  of  these  tissues.  Poulin  and  Krishnan 
(1996a)  algorithms  worked  well  only  for  fat:blood  PC  comparisons,  in  both  rats  and  humans. 

Performance  of  PCEB  versus  PCPB  Equations 

In  this  report,  PCEB  values  were  first  calculated  because  a  PBPK  modeler  would  use 
experimental  blood:air  numbers  when  these  data  were  available.  For  many  substances, 
however,  bloodiair  partitioning  experiments  have  not  been  performed  or  perhaps  the  results 
have  not  been  published.  Therefore,  PCPB  values  were  also  calculated  by  dividing  predicted 
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tissue:air  PCs  by  predicted  blood:air  PCs.  A  comparison  of  paired  (same  chemical,  species, 
tissue  and  source  of  algorithms)  PCEB  and  PCPB  values  are  shown  in  Table  10. 


Table  10:  Comparison  of  PCEb  and  PCPB  Tissue:Blood  Values  Calculated  from  the  Same 
Algorithms  for  each  Applicable  Tissue  and  Chemical 


Model  Type 

Number 

of 

Pairs 

Percent  of  Pairs  where 
Absolute  Percent  Error 
was  less  for  PCEB  values 
compared  with  PCPB  values 

Percentile  of 
Relative  Difference 

Min 

Median 

Max 

All 

103 

68 

2 

57 

153 

Empirical 

51 

55 

7 

35 

137 

Hybrid 

29 

66 

2 

58 

70 

Physiological 

23 

100 

74 

133 

153 

PCeb  values  tended  to  be  more  accurate  than  PCPB  values  overall.  Of  the  103  PCEB  and  PCPB 
paired  equations,  68%  had  a  smaller  absolute  percent  error  for  the  PCEB  value  compared  to  the 
PCPB  value.  Therefore,  tissuerblood  PCs  calculated  with  experimental  blood:air  data,  when 
available,  are  preferable  over  those  calculated  with  predicted  data  for  all  equation  types. 
However,  100%  of  PCEB  values  from  physiological  equations  had  smaller  absolute  percent 
errors  than  their  paired  PCPB  values.  This  indicates  that  predicted  blood:air  values  should  not 
be  used  to  substitute  experimental  blood:air  data  to  calculate  tissue:blood  PCs  from 
physiological  equations. 

In  half  of  all  pairs  (i.e.,  median  value),  PCEB  and  PCPB  values  differed  by  less  than  57%.  The 
range  of  difference  was  between  2  and  153%,  indicating  that  PCEB  and  PCPB  predictions  could 
be  very  similar  or  very  different.  Again,  it  is  clear  that  the  source  of  blood:air  value  used  for 
calculating  tissue.blood  PCs  from  physiological  equations  is  significant;  PCEB  and  PCPB  pairs 
were  at  least  (i.e.,  minimum  value)  74%  different  and  ranged  to  153%  difference. 


DISCUSSION 


Tissue:blood  and  bloodiair  partition  coefficients  for  each  tissue  compartment  are  essential  to 
PBPK  modeling  (Gargas,  1991).  Differences  in  tissue  and  blood  PCs  used  in  a  PBPK  model 
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can  influence  predictions  (Jang  etal.,  1999).  Overall,  solvent  PC  equations  predicting  rat  and 
human  blood:air  PCs  for  selected  JP-8  components  were  not  consistent  in  their  performance. 
Normal  biological  variation  in  blood:air  partition  coefficients  has  been  observed  at  20% 
(Gearhart  et  al.,  1993).  Of  all  predictions  (b!ood:air  and  direct  or  PCEb  tissue:blood),  just  30% 
were  within  ±20%  of  experimental  values.  It  is  important  to  note  that  all  types  of  equations 
(empirical,  physiological  or  hybrid)  performed  equally  well;  roughly  30%  of  calculations  from 
each  type  met  the  criterion  (±20%  of  experimental  value). 

The  complexity  of  the  tissue  in  question  appears  to  have  some  effect  on  which  type  of  equation 
performs  well  in  that  tissue.  Predictions  for  human  blood:air  PCs  were  more  tightly  clustered 
around  the  experimental  values  when  contrasted  with  rat  blood:air  PCs  (Figures  2  and  1 , 
respectively);  most  likely  the  higher  predictability  in  human  blood  is  due  to  decreased  capacity 
for  reversible  protein  binding  compared  to  rat  blood  (Poulin  and  Krishnan,  1996b).  Human 
blood:air  PCs  were  well  predicted  by  two  equations  (Gargas  et  al.  (1989),  Equation  23,  |20%|; 
Meulenberg  and  Vijverberg  (2000),  Equation  22,  |27%|)  whereas  rat  blood:air  PCs  were  well 
predicted  by  a  single  equation  (Gargas  etal.  (1988),  Equation  6,  |26%|  including  alkanes,  |14%| 
excluding  alkanes).  These  equations  were  either  empirical  or  hybrid;  all  were  derived  from 
experimental  data  for  low  molecular  weight  volatiles.  The  majority  of  chemicals  in  this  project 
were  chemically  similar  (i.e.,  relatively  small  volatiles),  making  an  empirical  equation  more  likely 
to  predict  the  experimental  data  well  (Payne  and  Kenny,  2002). 

Conversely,  fat  partitioning  is  relatively  simple  and  does  not  involve  significant  binding. 

Fat:blood  PCs  are  well  modeled  in  both  rats  and  humans  by  multiple  equations  (Table  9).  The 
best  of  these  were  physiological  equations  (Poulin  and  Krishnan  (1996a),  Equations  8  and 
29,|18%|  and  |14%|,  respectively).  The  success  of  these  algorithms  suggests  that  fat 
partitioning  depends  largely  on  the  composition  of  the  tissue  (fractions  of  neutral  lipids, 
phospholipids  and  water)  and  not  on  additional  factors  that  are  compensated  through  regression 
in  empirical  and  hybrid  equations.  Similarly,  Poulin  and  Krishnan’s  (1995b)  physiological 
equation  performed  reasonably  well  for  rat  muscle:blood  PCs  (Equation  50,  |16%|).  Fat  has  the 
highest  neutral  lipid  and  lowest  water  content  of  the  tissues.  Muscle  has  the  lowest  neutral  lipid 
and  highest  water  content  (Poulin  and  Krishnan,  1996a).  Physiological  equations  worked  best 
fui  the  two  most  differently  composed  tissues,  but  not  for  the  tissues  that  are  a  heterogeneous 
mixture  of  lipids  and  water. 
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Tissue:air  PCs,  paired  with  a  blood:air  PC  for  the  same  chemical,  tissue  and  species,  are  used 
to  calculate  a  tissue:blood  PC  using  Equation  1  (Fiserova-Bergerova,  1983;  Poulin  and 
Krishnan,  1995a).  Using  experimental  blood:air  PCs  (PCEB)  for  this  calculation  is  preferable 
since  68%  of  PCEB  values  generated  were  closer  to  experimental  data  values  (smaller  absolute 
percent  errors)  than  their  counterparts  that  were  calculated  from  predicted  bloodrair  PC  values 
(PCPB)  (Table  10).  Unfortunately,  experimental  bloodrair  values  are  not  always  available  in  the 
literature.  Although  vial  equilibration  studies  for  volatile  chemicals  are  not  difficult  or  expensive, 
nonvolatile  chemicals  require  expensive  radiolabeled  compounds  or  animal  intensive  testing  to 
determine  PCs  experimentally  (Murphy  et  al.,  1995;  Poulin  and  Krishnan,  1995a).  In  these 
cases,  PCPB  values  may  be  calculated  using  physiological  equations,  but  caution  should  be 
used  in  interpreting  the  results  since  100%  of  the  physiological  PCPB  values  were  farther  from 
experimental  data  as  compared  to  their  experimentally-based  PCEB  counterparts. 

Due  to  lack  of  solvent  PCs  for  nonane  and  decane,  only  one  comparison  was  possible  to 
determine  if  tissuerblood  PC  algorithms  were  suitable  for  alkanes.  Average  rat  bloodrair  PCs 
were  adversely  affected  by  the  inclusion  of  the  alkane,  decane,  with  the  exception  of  Gargas  et 
at.  (1989)  (Equation  5,  |44%  |  with  alkanes,  |51%|  without  alkanes).  Another  Gargas  et  al. 
equation  (1988,  Equation  6)  performed  better  (|26%|  including  alkanes,  |14%|  excluding 
alkanes).  In  general,  empirical  and  hybrid  equations  will  predict  better  for  chemicals  similar  to 
those  for  which  the  equation  was  derived  (Payne  and  Kenny,  2002)  and  are  of  limited  predictive 
value  beyond  these  chemical  families  because  they  do  not  incorporate  mechanistic  terms  for 
partitioning  (Poulin  and  Krishnan,  1996c).  The  results  of  this  with  alkane  vs.  no  alkane 
comparison  raise  doubts  as  to  whether  the  solvent  PC  based  algorithms  evaluated  in  this 
project  are  appropriate  for  the  task  of  calculating  tissue  and  blood  PCs  for  JP-8  components, 
which  are  largely  alkanes  (McDougal  et  al.,  2000). 

A  major  limitation  encountered  in  this  project  was  the  lack  of  experimental  data  for  the  selected 
JP-8  constituents.  This  deficit  included  solvent  PCs,  rat  tissue  and  blood  PCs  and  especially 
human  tissue  and  blood  PCs  for  these  chemicals  of  concern.  As  a  result,  predictive  algorithms 
for  only  four  rat  and  two  human  tissues  (blood  included)  were  able  to  be  adequately  evaluated 
in  Table  9.  Unfortunately,  this  lack  of  experimental  data  is  the  reality  for  scientists  wishing  to 
model  petroleum  based  mixtures.  The  most  prevalent  compound  (undecane,  Table  1, 

McDougal  et  al.,  2000)  in  JP-8  is  one  for  which  there  is  no  published  tissue  or  blood  partition 
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information;  however,  these  PC  experiments  in  rat  tissues  are  in  progress  at  the  University  of 
Georgia  (Fisher,  personal  communication,  2003). 

Although  the  lack  of  experimental  tissue  and  blood  PCs  for  these  chemicals  is  still  a  concern, 
additional  non-solvent  PC  based  algorithms  are  available  that  may  result  in  better  predictions  of 
experimental  tissue  PC  data.  These  algorithms,  however,  were  outside  the  scope  of  this 
project.  Alternative  algorithms  use  chemical  solubilities  (Paterson  and  Mackay,  1989;  Poulin 
and  Krishnan,  1996c)  or  other  physical  properties  such  as  geometric  volume  (Beliveau  and 
Krishnan,  2000). 


CONCLUSIONS 

Overall,  the  ability  of  all  three  types  of  equations  (Tables  2  through  5)  to  predict  blood:air  and 
tissue:blood  partition  coefficients,  both  necessary  to  PBPK  modeling,  was  limited.  Use  of 
calculated  PCs  from  solvent  partitioning  equations  should  be  carefully  considered.  For  the 
purpose  of  this  project,  use  of  experimental  PCs  (Tables  A-2  and  A-3)  in  the  JP-8  mixture  model 
is  recommended  whenever  available.  The  algorithms  shown  in  Table  9  may  be  used  when 
experimental  data  are  not  available.  For  all  applications,  if  calculated  PC  values  must  be  used, 
a  comparison  of  experimental  and  predicted  PCs  for  chemically  similar  compounds  is  advisable. 
This  allows  the  PBPK  modeler  to  account  for  the  variability  and  uncertainty  inherent  in  the 
partition  coefficient  used. 
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Table  A-1.  Solvent  Partition  Coefficients  for  Select  List  of  JP-8  Components 
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NA  -  not  applicable  -  The  temperature  dependence  of  log  P  is  very  small  (Leo,  personal  communication,  2003); 
however,  value  is  same  as  or  similar  to  25°C  value  in  most  of  the  instances  when  used  in  an  average. 


Table  A-2.  Published  Rat  Tissue  and  Blood  Partition  Coefficients  for  Select  JP-8  Components 
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Table  A-3.  Published  Human  Tissue  and  Blood  Partition  Coefficients  for  Select  JP-8  Components 
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Table  B-2.  Results  of  Algorithms  for  Human  Blood:Air  and  Tissue.Blood  (Direct,  PCEB  and  PCPB)  PC  Values 
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